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Abstract. The paper is a survey of some author's results re- 
lated with the Maslov-Shvedov method of complex germ and with 
quantum field theory. The main idea is that many results of the 
method of complex germ and of perturbative quantum field theory 
can be made more simple and natural if instead of the algebra of 
(pseudo)differential operators one uses the Weyl algebra (operators 
with Weyl symbols) with the Moyal ^-product. 
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Introduction 

This paper is a survey of some author's results related with the 
Maslov-Shvedov method of complex germ [2] and with perturbative 
quantum field theory. These results are, shortly, the following. 

Firstly, the results from the theory of quantum mechanical Schrod- 
inger equation. The Maslov-Shvedov method makes it possible to 
give a simple exposition of the method of canonical operator on a La- 
grangian manifold with complex germ [2-5] of asymptotic solution of 
the Cauchy problem for the Schrodinger equation. It turns out that 
many Maslov-Shvedov's formulas become more simple and get a nat- 
ural mathematical interpretation if in the Schrodinger equation one 
chooses the Weyl (symmetric) ordering of multiplication and differen- 
tiation operators. In particular, the transport equation describes the 
transport of half-forms along the classical trajectory, and the transport 
of Maslov-Shvedov wave packets is given by the action of the operator 
of the Weil representation of the metaplectic group corresponding to 
the tangent symplectic transformation to the Hamiltonian flow. Be- 
sides that, one obtains a simple and natural definition of the Maslov 
index modulo 4, related with the complex germ method at a point, 
which seemingly did not appear in the literature. A closed exposition 
of these results is given in §1. 

Secondly, the results from perturbative quantum field theory. Here 
introducing the infinite dimensional analog of the Weyl algebra allows 
one not only to interpret the quantum field theory Maslov-Shvedov 
complex germ obtained in [2] by complicated computations (or, more 
precisely, to obtain a result close to that result of [2]), but, generally, 
it allows one to give an exposition of the main results of perturba- 
tive quantum field theory not using subtraction of infinities from the 
Hamiltonian of a free field and normal ordering of operators. The Weyl 
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algebra plays the role of an extended algebra of operators in the Fock 
space, where the *-product corresponds to the composition of opera- 
tors, while the usual commutative product of functions corresponds to 
normally ordered product of operators. This result seems very impor- 
tant for understanding free and perturbative quantum field theory. A 
closed exposition of these results on the simplest example of the 
model of quantum field theory in four-dimensional space-time is con- 
tained in §2. 

The author is grateful to V. V. Dolotin, Yu. A. Neretin, and I. V. Tyu- 
tin for helpful discussions. 

1. Complex germ and the Weyl algebra in quantum 

mechanics 

1.1. The Schrodinger equation. The Schrodinger equation reads 

dtp d d 

(1) - ih— + H(t, -i/i— , - ih Q^-, Qi, ■ ■ ■ , Qn)*P = 0. 

Here ip(t, q±, . . . , q n ) is an unknown complex valued function (the wave 
function of a quantum mechanical system), H(t,p,q) is the Hamilton- 
ian of the corresponding system of classical mechanics, q = (qi, . . ., 
Qn), P = (pi, ■ ■ ■, Pn)- This equation is written in such a way that after 
substitution, instead of ip, of the quasiclassical asymptotics 

i S(t,q 1 ,...,q n ) 

(2) ij = a(t, qi ,...,q n )e l * 

(a and S are real functions varying very slowly when compared with the 
number h), in the principal approximation as h — > (the quasiclassical 
limit) one would obtain the Hamilton- J 'acobi equation 

OS TTf dS dS , n 

(3) - + /7(t,— , gi ,...,g„) = 

for the function S. 

The form of the Schrodinger equation (1) has some ambiguity in the 
general case. It is related with the fact that the operators of multi- 
plication by a function of qi and differentiation with respect to qi, in 
general, do not commute. Hence one should, in general, choose an or- 
dering of these operators in the quantum Hamiltonian H. In the case 
of a standard mechanical system without constraints, in which 

(4) H(t,p,q) = J2l^ + U(t,q), 

there is no such ambiguity. However, more general and deep consider- 
ations require to overcome this ambiguity. The usual way is to put the 
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operators of differentiation with respect to qi to the right of operators 
of multiplication by a function. 

1.2. Asymptotic Cauchy problem. Let us pose the asymptotic Cauchy 
problem for the Schrodinger equation (1): take oscillating initial data 



and let us look for an oscillating function ip(t, q) of the form (2), which 
turns into ipQ for t — and satisfies equation (1) up to o(h). To that 
end, two equations should hold: the Hamilton-Jacobi equation (3) and 
the transport equation, obtained by equating coefficients before h in 
the Schrodinger equation, into which the quasiclassical solution (2) is 
substituted. It is not difficult to see that, in the case when all the 
operators are put to the right of the operators q i: the transport 
equation reads 



Thus, our asymptotic Cauchy problem has been reduced to the Cauchy 
problem for the system of equations (3,6). It is well known that the 
Cauchy problem for the Hamilton-Jacobi equation amounts to inte- 
gration of the system of ordinary characteristic Hamilton differential 
equations 



(See, for example, [6], Ch. 4, or [1], Ch. 2.) Assume that this problem 
is already solved. How can one find the function a(t, q)l The transport 
equation is an ordinary differential equation for the function a, giving 
its behavior along the trajectories of the system of ordinary differential 
equations (7), where Pi = §f:- Considering particular cases (4) from 
standard quantum mechanics gives a solution of the form 



(5) 



MQi,---,Qn) = ao(q)e iS °W h , 




(7) 



dqi _ dH_ dpi, 
dt dpi ' dt 



OH 

dqi ' 



(8) 
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where (pi(t) , qi(t)) is a characteristic, i. e., a solution of the Hamilton 
equations. Indeed, let us differentiate ^jgy with respect to time: 

d %(t,g(0)) _ dHp. rj dp k ^ dqi 



(9) 



dt %(0) %(0) ^ mPk d qj {0) ^ mi %(0) 



_ if "9 2 ^ dqi 9gi 
2^ ^wp* /)„, /)„, a„ . + 2s n p 



dq k d qi dq 3 (0) ^ pm d qj (0Y 
whence 

da(t,q(t)) da ^-^ 9a 

= -2° {2^ H ^dq-dq- k +2s H ^)- 
This equation differs from the transport equation (6) by the term 

which is zero in the standard case (4). It turns out that here the 
point is the ordering of non-commuting operators —ih-£p and Oj in the 
Hamiltonian. If we choose a "right" ordering (recall that above we have 
arbitrarily put —ih-J^- to the right of qj), then the difference between 
equations (10) and (6) will disappear. The choice of a "right" ordering 
is the subject of the following two subsections. 

1.3. The Weil representation. As a simplest example consider the 
Hamiltonian H = pq with n — 1. For it, equation (10) requires the 
choice 

/-, 1 \ tt .. d ih ih ( d d 

(11) H = _^___ = __^._ + _. s 

instead of —ihq-j^ chosen above. That is, operators —ih-^p and qi 
should belong to the Hamiltonian symmetrically, without a prescription 
what stands to the right and what stands to the left. For formalization 
of these requirements, we need some information on the symplectic 
group. 

Consider the general quantum quadratic Hamiltonians 



>, n s fr v- 1 ih, f d d \ h 2 

12) H=^ -a jkmk - -h 3k \q 3 — + - -c jk 



_d__d_ 

dqj dq k 
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Here ajk and Cjk are symmetric real matrices, bjk is an arbitrary real 
matrix. In the Cauchy problem let us put the initial condition 

(13) ij = ij ZjP0 = ex P % - Q z am + J2 p o q ^j ■ 

Here Z = (Zij) is a symmetric complex matrix, p is a real vector. 
Assume that the matrix Z has positive definite imaginary part, then 
function (13) rapidly decreases at infinity (Gaussian wave packet). It 
turns out that, as it is not difficult to check by a direct computation, 
formula (8) gives in this case not only asymptotic but exact solution 
of the Cauchy problem: 

(U) ib(t a) = ^( az + b )( cz + d )- 1 >(( cz + d ) t )- 1 p° c ~^pZ(cz+d)-ic P0 
V ; ^ ' v/det(CZ + D) 

Here the sign T denotes transposing; A, B, C, D are matrices which can 
be found in the following way. The characteristics equations read 

dqi/dt = H Pi = ^2 bjiVj + c ^'' 

(15) j 3 

dpi/dt = -H qi = - ^ aijqj - Y bijPj 

3 3 

with the initial conditions 

(16) ft(0) = X>«?i(°)+P&- 

These are linear equations, hence, the evolution operator at the time t 
is a linear operator 



(17) expt 



-b -a \ _ f A B 
c b T -( C D 



In addition, this operator preserves the Poisson bracket of any two func- 
tions, as any evolution operator of the canonical Hamilton equations, 
i. e., it preserves the bivector field 

(18) " = £*- A ai' 

In the language of matrices this condition means that 

A B\ ( E\ ( A T C T \ _ ( E 
C D ) [ -E J ( B T D T J ~ ( -E 



(19) 



where E is the unit matrix. Such linear operators are called symplec- 
tic, they form the symplectic group Sp(2n, R). The set of symmetric 
complex matrices Z with positive definite imaginary part is called the 
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Siegel upper half-plane (cf. [7]); denote it by SQ. The group Sp(2n, R) 
acts on the half-plane SQ by the formula 

(20) Z i-> (AZ + B)(C Z + D)- 1 . 

Equivalently, this action can be defined by means of evolution at the 
time t of the matrix Riccati equation 

(21) Z + ZcZ + bZ + Zb T + a = 0, 

which is obtained by differentiating the action (20) with respect to t. 

Formulas (13), (14) define the action of the group Sp(2n, R) on the 
set of Gaussian wave packets. But this action is two-valued: to make 
it single-valued, one must choose one of the two continuous branches 
of the square root from det(CZ + D) ^ 0, Z e SQ. Hence a two-fold 
covering of the group Sp(2n, R) has a single- valued action on the set of 
Gaussian wave packets. This covering is called the metaplectic group; 
denote it by Mp(2n,R). 

It turns out that the action of the group Mp(2n, R) on the set of 
Gaussian wave packets is uniquely extended by continuousness to the 
action on the Schwartz space S = S(M. n ) of complex valued smooth 
functions ip(qi, . . . , q n ) rapidly decreasing at infinity, and also to unitary 
action on the space L 2 (R n ) of square integrable functions and to the 
action on the dual to S space S' = S'iW 1 ) of tempered distributions. 
Also the product of matrices yields the composition of operators. This 
representation of the group Mp(2n, R) is called the Weil representation, 
cf. [8,9]. 

The Weil representation is uniquely, up to a constant factor, char- 
acterized by the following property. Conjugation by an operator U 

corresponding to the matrix ^ ^ ^ ) ' P reserves ^ e 2n-dimensional 
vector space of operators with the basis 

, ., d d s 

(22) (q 1 ,...,q n ,ih—,...,th- 



dqi dq n ' 

( A B 

and acts on this space by the matrix I „ n I . This is obtained 



x C D 

by exponentiating from the fact that the commutator with the Hamil- 
tonian (12) also preserves this space, and acts on it by the matrix 

, y J , up to the factor ih. This property implies the unique- 



c b 1 

ness of the Weil representation as follows. If U' is another operator 
with the same property, then the operator U'U^ 1 commutes with the 
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operators qj and ih-£p, and hence it is multiplication by a constant, as 
it is not difficult to show. 

In particular, the matrix ^ ^ ^ ^ acts by multiplication by the 

function exp J2j ;,k BjkXjXk) the matrix ^ ^ {A^)' 1 ) acts ^ y 
composition of a linear change of coordinates given by the matrix A, 
and multiplication by V det A. Finally, the matrix ^ ^ ^ ^ acts 
(up to a constant factor) by the Fourier transform: 

(23) (F h i,)(q) = ^= / e^^^dy. 

This transform exchanges the operators g,,- and (up to sign); the 

square of this transform is the change of variables q — > —q. In quan- 
tum mechanics the Fourier transform of the wave function is called its 
momentum representation. 

The above matrices generate the group Sp(2n, R), which gives a proof 
of existence of an action of the group Sp(2n, R), defined up to a factor, 
with the above described commutation relations with the operators 
(22). 

Note also that exponentiating of operators (22) yields an action of 
the so called Heisenberg group H n on the space of functions. Namely, 
the operator a\qi + . . . + a n q n acts by multiplication by the function 
e l ^ a i q il h , and the operator ih^bi-J^ acts by the change of variables 
q% — ► q% + bi. The multiplication in the group H n is defined with the 
help of the formula 

(24) e a e S = e h+h ■ e c/2 , ab-ba = c 

(here c is a number). This action is compatible with the action of the 
group Mp(2n, R) in an obvious sense, so that the space of functions 
has an action of the semidirect product of the groups Mp(2n, R) and 
H n . This group is a central extension of the affine symplectic group 
(i. e., the semidirect product of the group Sp(2n, R) and the group of 
parallel translations in the space R 2n ) with the help of the circle. Let 
us denote the action of an element g of any of these groups (possibly 
defined up to a factor) on the space of functions by the symbol 

U = P (g). 

1.4. The Weyl calculus. This is a way to assign to a function (p(pi, 
. . ., p n , qi, . . ., q n ) an operator ip = ip{p, q) on the space of functions 
ip(qi, . . . , q n ), this correspondence being in accordance with the action 
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of the affine symplectic group, i. e., for any element g of this group we 
have 

(25) p(g)tp(p, q)p(g)~ 1 = gp(p, q), 

where (g<p)(p,q) = y^fi 1 " 1 ^, <?))• The correspondence ip — > <p possesses 
also the following properties: 

a ) ((I2 a iqi + b iPi) T= (E^ -^^-) fc ; 

b) a real function ip corresponds to a symmetric operator <p, i. e. 
such that 



ip 1 ■ (pip 2 dq = J (ftp! ■ ip 2 dq 



for any rapidly decreasing functions ipi, ip2', 

c) appropriate continuousness properties, into which we shall not go, 
see, for example, Hormander's book [10]. 

Using these properties one can define the operator ip> for a rather wide 
class of functions <p. First of all, for polynomials ip(p, q) the operator 
(p is defined uniquely from property (a). For example, 

(26) (pqf= -((p + qf - p 2 - q 2 y = -ih [q— + - 

Similarly, in algebra a way is known to express each polynomial of 
2n variables through powers of linear forms. For any homogeneous 
polynomial of degree k of 2n variables there exists a unique symmetric 
/c-linear form of 2n variables (the polarization of the polynomial), which 
gives this polynomial for coinciding arguments. A polylinear form is 
a tensor of rank k, i. e. an element of non-commutative algebra of 
2n generators. Substituting instead of these generators the operators 
qi and —ih-^-, we obtain the required operator. This operation is 
GL(2n, M) -invariant. 
Further, we have 

(27) (^xp^C^aiPi + b i q i ) S j^= exp (^2^-^- + ^>%<hj ■ 

Since many functions ip can be expressed as superposition of exponents 
of linear forms using the inverse Fourier transform: 

(28) <p(p, q) = J (F h <p)(a, &) e *E«*«+fc*> dadb 

(see (23)), we obtain a way to find the required operator for a large 
class of functions. Let us give the explicit formula for the operator <p: 

(29) (W)(q) = J J <p(p, (q + y)/^^) d pdy. 
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Further, it is not difficult to compute the formula for the * -multiplication 
(fii*<f2, i- e., for the function corresponding to composition of operators 
(fii(f>2, so that 

(Vl * ^2)"= <£lV?2- 

For that one should compute the product of two operators of kind (27) 
by formula (24), and use the inverse Fourier transform (28). Let us 
give the answer. Denote 



(30) { Vl , V2 } = J2 



where = g« for 1 < i < n and y t = Pi- n for n + 1 < i < 2n, and 
u %3 = Sij^n - 5i_ n j. Then 

(31) (<pi * <P2)(Vi) = exp (-yX!^ /.^. ) PMwizi 

— i/l 

This product is usually called the Moyal product; it is not difficult to 
check directly that it is associative. 

Finally, the same formula (24) implies that the transport equation 
has the right form (10) for the Hamiltonian being exponent of a linear 
form, and hence by linearity for any Hamiltonian. Below the Hamilton- 
ian H and other quantum observables will be understood in the sense 
of Weyl calculus. 

1.5. Method of complex germ at a point. Following Maslov and 
Shvedov [2], let us look for asymptotic solutions of the Schrodinger 
equation (1) in the form of wave packets 

(32) ij(t, q) = f (t, ^=^) e K£i*<«><*-«&<*>H*<*>) 

for some functions q l (t), p l (t), S(t), f(t,x). Let us call them Maslov- 
Shvedov wave packets. An example of such wave packet is the Gaussian 
wave packet (13). We will find equations on these functions which will 
imply that the wave packet (32) satisfies the Schrodinger equation up 
to o(h). 

To this end, note that: 

1) action of the operator fa — q l on the function ip amounts to action 
of the operator \fhxi on the function /; 

2) action of the operator pi —p % on the function ip amounts to action 
of the operator —iy/h-^ on the function /; 
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3) the Schro dinger equation (1) for the function ip amounts, up to 
o(h), to the equation 

(5>o?o - s)f + VhJ2(-p>* - <-^)f + ^% 

d 

= H(p ,q )f + Vh^2(H q .Xi - iH p .—)f 

hf^ f _8_ _d_\_rr JM 

+ 2 [2^ H ^ x i %ti q ^ [Xi + dx _ XiJ H PiPj ^ ^ J } 

(all derivatives of the Hamiltonian are taken at the point (p , q )). This 
equation will be satisfied provided the following system of equations 
holds: 

(33) 

^JPoQo -S = H(p ,q ), 

<?o = H pv 
Po = ~H qi , 

■?L--(\^rr -m ( — JL\_h ——\f 

% 8t ~ 2 \ l^ H ^ XiX i lH ™ \ Xi dxj + dxj Xi ) ^dxidxj) /- 

These equations mean that (po(t) , q (t)) is a classical trajectory, S(t) 
is the action along this trajectory, and the function / satisfies the 
Schrodinger equation with the quadratic Hamiltonian depending on 
time, with h — 1. 

This latter equation means the following. The quadratic part of 
the Hamiltonian at each point of the trajectory gives an infinitesimal 
symplectic transformation: 

dxi/dt = H qjPi Xj + Y HpiPjVj, 

(34) / 

dyi/dt = ~2_^ H qiqjXj - 2^ H qiPjVj . 

3 3 

The composition of all these transformations at the time t gives a 
metaplectic transformation 

(35) *® ),Vdet(C7(0Z + £>(«))), 

whose action on the function f(0,x) under the Weil representation 
gives the function f(t,x). 
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In particular, if we look for the function / in the form of a Gaussian 
function 

(36) /(( - x) = Tsmrn + w> exp G - z>i(t)XiX ') • 

then for the function 

(37) Z(t) = (A(t)Z(0) + B(t))(C(t)Z(0) + D(t))- 1 
we get a matrix Riccati equation 

(38) Z + ZH PP Z + H pq Z + ZH* + H qq = 

of type (21). The matrix is called the complex germ. 

It is rather interesting to express these equations in terms of the 
functions q(t), q'(t) and the Lagrange function F(t, q, q'), i. e., to rewrite 
them in the language of the variational principle. Then equation (34) 
turns into the Jacobi equation in the theory of second variation, and 
equation (38) turns into the corresponding matrix Riccati equation, see 
Gelfand-Fomin's book [11], cf. [12]. But in the variational calculus the 
matrix Z is real; this case will be considered below. 

1.6. Method of canonical operator. In conclusion we shall briefly 
discuss the powerful method of canonical operator, due to V. P. Maslov. 
This method allows one, for example, to write out the asymptotic solu- 
tion of the Cauchy problem for the Schrodinger equation (see 1.2). For- 
mulas (2), (8), (3) yield this solution for t sufficiently small, when dif- 
ferent characteristics do not intersect each other and det ( j^^ j ^ 0. 

Method of canonical operator shows what happens with the solution 
after passing through focal points, where characteristics intersect each 
other and the determinant vanishes. To this end, let us represent the 
solution (2) as a superposition of wave packets (32) satisfying equations 
(33): 

(39) Mq) = f e^W+^X^))); ^ ^J^) ^, 
a = (a 1 , . . . , a n ). 

Let us first consider the case when the n-dimensional submanifold 
(po(a), qo(a)) of the phase space (p,q) diffeomorphically projects onto 
the g-plane. Let us develop the expression under the exponent into the 
Taylor series in a vicinity of the point a for which qo(a ) = q, and let 
us make change of variables 

ao — a 

~7T = V ' 
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We obtain 
(40) 



S(a)+p (a)(q - q (a)) = S(a ) + VhJ^ [ + 



Hence ^(g) will not decrease more rapidly than any power of ft, only in 
the case when 



Below we will assume that this equality holds identically. It implies 
that 

(42) E (MM_MM) =0 

for all In other words, the symplectic differential 2-form 

(43) lo = dpi A d^j 

vanishes on the submanifold (po(a), qo(ct)) of the phase space. Maslov 
called such submanifolds Lagrangian. Conversely, any Lagrangian sub- 
manifold diffeomorphically projecting onto the g-plane, is the graph of 
the differential of some function S(q): 

(99 

(44) Pi(?)=«-, l<i<n. 
Equality (41) also implies that 

iS(ar\) 

J e -^ TpQ ' lx f(a ,x)dx + O(Vh) 



e h 



| detQ| 
e l5(<?)/fe a(g) +0(v / ft), 



14 



A. V. STOYANOVSKY 



Let us now assume that the Lagrangian submanifold has been trans- 
formed by the Hamiltonian flow (7) on the phase space at the time t. 
What then happens with the functions S(q) = S(ao) and a(q)7 

Recall that on each trajectory the tangent metaplectic transforma- 
tion (35) arises. Denote it by g. Then in the formula (45) the following 
changes will occur: 

1) po, qo are transformed by the flow; 

2) S(ao) — > 5(0:0) — 5'(ao)+ a ction along the trajectory; 

3) P -> AP + BQ, Q -> CP + DQ; 

4) f(a ,x) -> p(g)f(a ,x). 

How is a(g) transformed? To answer this question let us study what 
is p(g)ei x PQ x , as promised at the end of 1.5. 

Denote Z = PQ^ 1 . Now Z is real. Moreover, ipz(x) = e% xTZx does 
not belong now to the Schwartz space and to L 2 . It is the unique, up 
to proportionality, distribution solution of the system of equations 

(46) (i^+£^Wz = 0, 

The transformation g takes these equations to the equations 

(47) ^ (vijXj + w ij i -^r^j ^ = 0, 1 < % < n, 

which, for det(CZ + D) ^ 0, are equivalent to the equations on the 
Gaussian function iP(az+b)(cz+d)- 1 ■ Moreover, as it is shown by taking 
the limit in equalities (13), (14) as ImZ — > 0, po — (ImZ is the 
imaginary part of the matrix Z), for det(CZ + D) 7^ we have 

^ = ..„l n ; + o Vdet( c Z+g) »'^'^'- = 

( 48 ) p mfc/2 



VI det(CZ + D; 



=V ; (AZ+B)(CZ+Z))- 



for some integer k called the Maslov index. This index has a purely 
algebraic definition, see, for example, Hormander's book [10], §21.6. 
As far as we know, the above simple definition of the Maslov index did 
not appear in the literature. 

In the general case (when det(CZ + D) can equal 0) it is easy to see 
that the system of equations (47) forms a basis of a (real) Lagrangian 
subspace L in the 2n-dimensional vector space with the basis 

, A . d d 

(49) [xi,...,x n ,i— ), 

OXi ox n 
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i. e., an n- dimensional subspace on which the symplectic form, given 
by the commutator of operators, vanishes. It is not difficult to see 
that the solution ip(x) = of this system is, in general case, up 

to a constant factor, the product of a function of type ipz of part of 
the variables, for a real Z, and the delta function of the remaining 
variables (after an appropriate linear change of coordinates x). The 
most degenerate case is the system of equations x^ — 0, 1 < i < n, 
whose solution is the delta function S(x). 

Thus, we have described the Mp(2n, M) -orbit of functions ipz for real 
Z, or, which is the same, the orbit of the function 1 in the projectiviza- 
tion PS' of the space S' of distributions. This orbit is isomorphic to 
the variety of real Lagrangian subspaces of the 2n-dimensional symplec- 
tic vector space. This variety is called the Lagrangian Grassmannian; 
denote it by A n . The embedding A n — > PS' induces a complex line 
bundle /i on the Grassmannian A n , whose fiber at the point L e A n is 
the line O0£. This bundle has an action of the group Mp(2n, M). Triv- 
ializations and transition functions of this bundle can be obtained from 
formula (48). Let us call this bundle \i the Maslov bundle (Hormander 
[10], §21.6, uses another terminology and calls the bundle \x the tensor 
product of the half densities bundle and the Maslov bundle). 

Returning to formula (45), we see (assuming that f(ao,x) belongs 
to the Schwartz space with respect to x) that under the action of the 
Hamiltonian flow the function a(q) is multiplied by 



if after the transformation by the flow the Lagrangian manifold still 
diffeomorphically projects onto the g-plane (i. e., if det (CP+DQ) ^ 0). 
This result generalizes formula (8). 

Let us now return to the integral (39) and consider it in the case when 
the manifold (po(a), qo(a)) not necessarily diffeomorphically projects 
onto the g-plane. In this case let us present the function f(a,x) as 
sum of functions fi, each of which has the support with respect to the 
variable a diffeomorphically projecting onto some Lagrangian plane in 
the phase space. Then let us apply to these functions linear Hamil- 
tonian flows (15), giving metaplectic transformations gi, so that these 
Lagrangian planes get to the g-plane. After that let us apply formula 
(45). We obtain a wave function tpi(q). Finally, let us apply to these 



det Q\y/\ det(C PQ- 1 + D)\ 
e™ k ' 2 1 det (CP + DQ) \ 



^\o\et(CPQ- l + D)\ 



e 



ink/2 



(50) 



-iirk/2 
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wave functions the transformations p(gf 1 ), and let us take the sum of 
them. The obtained wave function ip(q) (which is, in general, a dis- 
tribution) is defined correctly up to 0{\fh) if the Maslov index of any 
closed curve on the Lagrangian manifold is divisible by 4. The role of 
function a(q) is played here by a section of the Maslov bundle on the 
Lagrangian manifold, induced from the bundle \i on the Lagrangian 
Grassmannian. 

Thus, the method of canonical operator assigns a (distribution) wave 
function ip(q), defined up to O(Vh), to a Lagrangian submanifold and 
to a section of the Maslov bundle on it. In the case when the La- 
grangian submanifold is the graph of differential of a function S(q), 
the wave function has the form (45). Under the evolution given by 
the Schrodinger equation, the corresponding Hamiltonian flow on the 
phase space transforms the Lagrangian manifold and the section of the 
Maslov bundle on it, and hence transforms the wave function i/j(q). 
This gives the global asymptotic solution of the Cauchy problem. 

The method of canonical operator has far generalizations. For exam- 
ple, if one integrates not along n-dimensional but along /c-dimensional 
(k < n) isotropic submanifold, then one obtains the method of canoni- 
cal operator on a Lagrangian manifold with complex germ. The method 
of complex germ also applies to approximate solution of linear and even 
non-linear partial differential equations, and not only to the Schrodinger 
equation. See Maslov's books [3-5]. 

2. The Weyl algebra in quantum field theory 

2.1. The Schrodinger equation. Consider a field theory action func- 
tional of the form 



dx n , 



(51) J= f F(x°, . . . , x n , u 1 , . . . , u m , u].o, . . . , u™„) dx° . . . 

JD 

where x° — t, the independent variables, u 1 , . . . ,u m are 

the dependent variables, u l xj = |^-, and integration goes over an (n + 
l)-dimensional surface D (the graph of the functions u l (x)) with the 
boundary 3D in the space K m + n+1 . The main simplest example for our 
considerations is the (p 4 model in four dimensions: 

(52) J = j ^ (u 2 t - Yl u l j ~ m ^ 2 j - \ 9U ^j dtdxldx2dx3 - 
The Schrodinger equation for the model (51) reads 

(53) ,h- = jH[t, x, „'(x), -ih-^ j * dx. 



MASLOV'S COMPLEX GERM AND THE WEYL-MOYAL ALGEBRA 17 



Here x = (xi, . . . ,x n ); 

\& is an unknown complex valued functional of the variable t and of 
functions w*(x), 1 < % < to; 

iif is the density of the Hamiltonian of the theory, which equals the 
Legendre transform of the Lagrangian F with respect to the variables 
u\\ denote the dual variables to u\ by p l ; 

H is the density of the quantum Hamiltonian, obtained from H by 
the substitution of the variational differentiation operator ~ih s J^ 

instead of p l . 

As in quantum mechanics, here the problem arises of ordering of the 
operators w l (x) and —ih g J^ in the quantum Hamiltonian. Let us not 

consider this problem now, all the more in the case of the ip A model 
there is no such ambiguity. 

The Schrodinger equation (53) is obviously relativistically non-invar- 
iant. But one can write out a relativistically invariant version of this 
equation, in which the surface t = const in the space-time is changed 
by an arbitrary space-like surface, and the functional \l/ depends on 
this surface and on functions u l (s) on it, where s = (si, . . ., s n ) are 
parameters on the surface. This relativistically invariant version is 
obtained in exactly the same way as the usual quantum mechanical 
Schrodinger equation, by the formal substitution into the generalized 
field theory Hamilton- Jacobi equation, see [6]. In physical literature 
a close equation is called the Tomonaga-Schwinger equation [13], and 
the problem of solving this equation is called quantization on space-like 
surfaces. 

One can give a rigorous mathematical sense to the Schrodinger equa- 
tion (53) and its relativistically invariant generalization. To this end, 
the conventional usual way is to consider weakly continuously differ- 
entiable sufficient number of times functionals on a nuclear space of 
functions u % (s), for example, on the Schwartz space. In this interpreta- 
tion, for example, the variational derivative s ^ x ^ is a distribution in x, 

the second variational derivative f r, ttv, » is a distribution in (x, x'), 

bu l (x)bu l (x'j v ' ' 

etc. 

However, it is not difficult to see that with such understanding the 
Schrodinger equation, say, for the tp A model, 
(54) 

ih^ = j {- Y J^ + 2(grad.(x)) 2 + ^u(,f + -.(x) 4 j * dx, 

does not have nonzero four times differentiable solutions. Indeed, con- 
sider the second derivative In the expression for this derivative 
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following from equation (54), we will have the term 

/ / ^" (y)2 * <ixrfy - 

which, as it is easy to see, has no sense (the second variational derivative 
cannot be restricted as a distribution to the diagonal x = x'). 

One can give physical arguments as well in favor of the statement 
that states cannot be functionals. Indeed, if it were so (as it was 
assumed in past, see, for example, [14,15], etc.), then the values of 
these functionals or related quantities, in principle, could be measured. 
On the other hand, it is known (see, for example, §1 of the book [16] 
by Berestetsky, Lifschitz, and Pitaevsky) that in relativistic quantum 
dynamics, quantum mechanical quantities like energy and momentum 
are theoretically non- measurable, and the only measurable quantities 
are the scattering sections. 

Besides that, one would like to have that in the case of free scalar 
field given by a quadratic Hamiltonian, the Schrodinger equation be 
solved exactly, similarly to the finite dimensional case. This implies 
that the space of states and the space of operators have an action of 
infinite dimensional symplectic group. Indeed, the evolution operators 
of classical field equations from one space-like surface to another are 
canonical transformations, preserving the field theory Poisson bracket. 
This follows from the generalized field theory canonical Hamilton equa- 
tions (see [6]). In the case of free field these operators are linear, i. e., 
symplectic. Hence the action of quantum Hamiltonians should admit 
a compatible action of a group of symplectic transformations of the 
space of functions (u l (s), p l (s)). 

The traditional action of an infinite dimensional symplectic group 
is the projective Segal-Shale-Weil-Berezin representation in the Fock 
space [14,15]. However, this infinite dimensional symplectic group does 
not suit for our purposes, as shown in the important paper [17]. In this 
paper it is shown that the evolution operators of the Klein-Gordon 
equation from one space-like surface to another, in general, do not 
belong to that version of infinite dimensional symplectic group which 
acts on the Fock space. This is also in accordance with the physical 
arguments above. The evolution operators belong to the group of con- 
tinuous symplectic transformations of the Schwartz space of functions 
(u l (s) , p % (s)) . It is this group that should act on the space of quantum 
Hamiltonians. 

To achieve this, it is natural, as in §1, to introduce, instead of the 
algebra of differential operators on the space of functions, the infinite 
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dimensional generalization of the Weyl algebra. Let us give a definition 
of this generalization. 



2.2. Infinite dimensional Weyl algebra. 



2.2.1. Definition of the Weyl algebra. The Weyl algebra is constructed 
starting from a symplectic vector space. Consider the symplectic Schwartz 
space of rapidly decreasing functions (V (s) , p 1 (s)) with the Poisson 
bracket 

(55) $2} " \ J {s^(s)Sp^s) " Wis) Ms) ) 

of two functionals $;(«*(•), p*(-)), I = 1,2. Let us write it in the form 

where y 1 = u l for 1 < i < m and y l = p l ~ m for m + 1 < % < 2m, 
and cj u = 8ij- m — 5j_ m j, as in 1.4. The Weyl algebra is defined as 
the algebra of infinitely differentiable functionals $(V(-), p l (-)) with 
respect to the Moyal *-product 
(57) 

($ 1 *$ 2 )(y i (-)) 



«*(•)=»*(•) 




This product is not everywhere defined: for example, u l (s) * p l (s) is 
undefined. We shall not go into details of the domain of multiplication, 
as well as into details of defining topology in the Weyl algebra. This 
should be the subject of a separate investigation. Note only that if all 
necessary integrals and series are defined and absolutely convergent, 
then the *-product is associative. This is a formal check similar to 
the finite dimensional case. In this paper we will be interested only in 
some concrete computations in the Weyl algebra. In algebraic quantum 
field theory [17,18,19] a somewhat different definition of Weyl algebra 
is adopted. 

Below we will see that the Weyl algebra allows one to construct 
a logically self-consistent theory of free quantum scalar field and to 
simplify drastically perturbative theory of interacting quantum fields. 
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2.2.2. The problem of states. Thus, operators in equations (53), (54), 
and others will be understood as elements of the Weyl algebra. And 
how will be understood states VP? They already cannot be function- 
al of u l (-), since the Weyl algebra does not act on them. In the case 
of a finite dimensional symplectic vector space, the Weyl algebra acts 
canonically on half-forms on a Lagrangian subspace. In coordinates 
qi,..., qN,Pi, • • • ,Piv half-forms look f(q 1 , . . ., q N ) (dqi . . . dq N ) l/2 . The 
fact that the Weyl algebra acts on half-forms, can be seen, for exam- 
ple, as follows. The operator f(pi?j)T of infinitesimal linear change of 
coordinates from the Weyl algebra acts as 



and this is the action on half-forms. 

What are half-forms on an infinite dimensional space of functions 
■u l (s)? Seemingly, one cannot say anything definite at this point. At 
least, half- forms cannot be constructed from finite dimensional spaces, 
analogously to the construction of measures on an infinite dimensional 
space. Author's attempts to construct half-forms failed to be success 
(see, for example, [20]). 

But actually, in order to obtain physically important quantities for 
free field, we need not states: it suffices to use only operators, as will 
be shown below. States are "non-observable neither physically nor 
mathematically". Hence we will consider, instead of equations (53), 
(54), the Heisenberg equation for an element u l (-),p l (-)) of the 
Weyl algebra: 



is the commutator in the Weyl algebra. The classical limits of equations 
(58) are the field theory Hamilton equations 




(58) 




and its relativistically invariant generalization, where 



(59) 



[$i, $ 2 ] = $1 * $2 - $ 2 * $1 



(60) 




equivalent to the Euler-Lagrange equations. 



2.3. Quantization of free scalar field. 
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2.3.1. Solution of the Heisenberg equation for free scalar field. Solution 
of equation (58) is given by the formal equality 

(61) $(f0 = U(t , h) * $(f ) * U(t , h)- 1 , 
where 

(62) U(to,h) = Texp f [ ^-H(t,x) dtdx, 

Jto J lh 

and Texp J means the ordered exponent (the multiplicative integral): 
(63) 

Texp J T(t)dt=l+ J F(t)dt+ J T(t) * T(t') dtdt' + . . . . 

<0 t <t<ti to<t'<t<tl 

(Cf. 2.4.1 below.) 
Let us first consider the free scalar field (g — 0). In this case 

if(f,x) = -(p(x) 2 + (gradw(x)) 2 + m 2 w(x) 2 ) 

is a quadratic expression not depending on t, hence we can omit the 
sign T before exponent. Due to the fact that the Hamiltonian 

(64) H = J H(x)d: 
is quadratic, we have 

(65) i-[tf ,$] = {$,#o}, 

therefore, ®(h; u(-),p(-)) is obtained from &(t ; u(-),p(-)) by the linear 
symplectic change of variables 

(66) (u(i ,x),p(t ,x) = u t (t 0: x)) -> x),p(ti,x) = ttt(ti,x)) 

given by the evolution operator of the canonical Hamilton equations, 
i. e., by the evolution operator of the Klein-Gordon equation 

(67) Du - m 2 u = -— + ^ ^-^ - m 2 « = 

from the Cauchy surface t — to to the Cauchy surface t — t 1 . Here 
it is rather essential that the evolution operator is a continuous linear 
invertible operator in the Schwartz space of functions (w(x), p(x)). 
Similar statement is true for the evolution of the functional $ between 
any two space-like Cauchy surfaces. (For non-quadratic Hamiltonians 
and non-linear classical evolution operators similar statement is not 
true.) 



x 
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Hence we can identify the Weyl algebras corresponding to different 
space-like surfaces, by means of the evolution operators of the Klein- 
Gordon equation. (Cf. [17].) In other words, we can consider the Weyl 
algebra Wq of the symplectic vector space of solutions w(t,x) of the 
Klein-Gordon equation on the whole space-time. The symplectic form 
on this vector space is given by taking the Cauchy data 

(68) u(t,x)->(u(s),p(s)) 

on any space-like surface x = x(s). (The quantity p(s) is proportional 
to the normal derivative of the function w(t,x) at the point s.) Below 
we will fix this identification of the Weyl algebras of various space-like 
surfaces. 

2.3.2. Green functions. Let us now consider the free scalar field with 
a source, i. e. put 

(69) H(t,x,u,p) = -(p(x) 2 + (grad M (x)) 2 ) + -y«(x) 2 +j(*,x)w(x), 

where j(t,x) is a smooth function with compact support (a source). 
Denote the corresponding formal element (62) of the Weyl algebra by 
Uj(to,ti), and the Hamiltonian by 

(70) H s (t) = H + J j(i, x)u(x) dx, 

to show dependence on the source. Then, if the support of the function 
x) is situated between the planes t = t min and t = t max , then the 
formal element 

(71) -^j(^o) — Uo(t mSuX ,to) * t/j(£ m i n , tmax) * C^o(^O) ^min) 

of the Weyl algebra does not depend on t min ,t mSuX . Besides that, we 
have 

(72) R^h) = U (t ,t 1 )*R i (t Q )*U (t ,t 1 )- 1 . 

Hence the element -Rj(t ) = Rj(t ; u(-) , p(-)) correctly defines an ele- 
ment of the Weyl algebra of any space-like surface under our identi- 
fication, i. e. an element Rj of the Weyl algebra W . This element 
equals 

oo 

(73) R J = R i (u(-,-))=Texp J J -^j(t,x)u(i,x) dtdx, 

— oo 

where 

(74) u(t,x) = exp (~~~j~ H oj * M ( x ) * ex P f " Jj~ H o 
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is understood as a functional on the space of solutions «(•,•) of the 
Klein-Gordon equation, i. e. as an element of the algebra W . (See 
2.4.1 below.) Let us emphasize that expression (74) is purely symbolic, 
since the element exp(£i —t )H /(ih) does not exist in the Weyl algebra, 
because already H * H does not exist. Let us call the element Rj the 
generating functional of operator Green functions of a free field. Let us 
also call the coefficients of the Taylor decomposition of the functional 
Rj with respect to j at the point j = 0, 



(75) (ih) 



n $ R i 



5j(ti,xi) . . . 5](t N ,^ N ) 



= Tu(t!, Xi) * ... * u(t N , Xjv), 



by the operator Green functions of a free field; here the symbol T de- 
notes *-product ordered by decreasing of the variables tj. The operator 
Green functions are distributions of (ii,Xi), . . ., (£tv, x tv) with values 
in W , symmetric with respect to permutations of indices. 

Let us now pass to the scalar Green functions. To this end, define a 
linear functional on the algebra Wo, called the vacuum average of an 
element $ from Wo and denoted by (<&) or (0|$|0), in the following 
way. The momentum representation 

(76) u(p , . . . ,p n ) = {27r) l +1)/2 J < ' X >' </(.r". ...,x n )dx 

of a solution u(t, x) of the Klein-Gordon equation, where t = x°, 
is a distribution supported on two sheets of the mass surface p = 
±a/p 2 + m 2 , where p = (pi, . . . ,p n ); in this paper we restrict our- 
selves by theories with nonzero mass, m > 0. Hence w(t,x) can be 
uniquely decomposed into the sum 

(77) u(t, x) = u+(t, x) + u-(t, x) 

of a positive frequency solution u + (t,x), whose Fourier transform is 
supported on the sheet po > 0, and a negative frequency solution 
u-(t, x), whose Fourier transform is supported on the sheet p < 0. 
We have 

[u+(fl,Xl),U + (*2,X 2 )] = [U-(fl,Xi),U_(i 2 ,X 2 )] = 0, 

[u4p),u + (p')] = -h5(p + p')5(p 2 -m 2 ), p <0, Po >0, 

where p 2 = pi — YTj=iP 2 j- Define ($) as the unique (not everywhere 
defined) functional with the following properties: 

(79) ($*u_(i,x)) = («+(*, x)*$) = 0, (1) = 1. 
Define the Green functions by the equality 

(80) Xi) . . . u(t N , Xjv)) = (Tufa, Xi) * ... * M^at, Xjy)), 
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and their generating functional by the equality 

(81) Z(j) = (Ri). 

A computation left to the reader (apply Fourier transform with respect 
to x; cf. the textbook [21] by Bogolyubov and Shirkov) shows that 
the two-point Green function turns out to be equal to the Feynman 
propagator 

(82) (u(t, X )u(t>, x')>~ = ih 2 S(P + /1 ■ . 

and the generating functional of the Green functions is given by the 
usual expression 

83 Z(j) = exp — / — — — dp. 

2h J p 2 — m 2 + ie 

2.3.3. The Fock space. Define the standard Fock space, linearly gener- 
ated by the vectors 

(84) |p ( i) , . . . , p (N) ) = u + (p ( i) )...u + (p (N) ) 1 0) 

(after integration over p^ with a function f(p(i), ■ ■ ., P(n))) f° r all N, 
. . ., p(jv) such that = m 2 , p (i) > 0. In other words, the Fock 
space is the direct sum over all iV of spaces of symmetric functions 
of iV variables puy On this space one introduces the structure of a 
Hilbert space, namely, the direct sum over all iV of the spaces L 2 of 
symmetric functions of iV variables p^) with respect to the natural 
Lorentz-invariant measure 

(85) 5(p 2 - m 2 )dp = 

on the mass surface p 2 = m 2 . 

One can formally assign an operator in the Fock space to an element 
$ of the Weyl algebra Wq, with the matrix elements 

(86) (0|u_(-p' (1) ) . . .U-(-p[ NI) ) * $*«+(p ( i)) . . .u + (p (N) )\0). 

But for many important operators responsible for local dynamics, for 
example, for the Hamiltonian $ = H , the expression (86) is undefined. 

Note two properties of this correspondence, which it is not difficult 
to check. 

1) ^-product of functionals corresponds to composition of operators, 
so that this correspondence is a (not everywhere defined) homomor- 
phism of the Weyl algebra W to the algebra of operators in the Fock 
space. 
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2) Complex conjugation of functionals goes to the Hermitian conju- 
gation of operators in the Hilbert space. In particular, the operator 
u+ip) is Hermitian conjugate to the operator u + (p) = U-(—p). 

2.4. Quantization of interacting fields. We start with the formal 
decomposition of a solution of the Heisenberg equation (58) in the 
model into the perturbation series with respect to the coupling con- 
stant g. For that, recall the perturbation theory of linear differential 
equations. 

2.4.1. Perturbation theory of linear differential equations. Consider the 
equation 

(87) ^ = A(t)v(t) + B(t)v(t), 

where A(t) is a linear operator applied to a vector v(t), and B{t) is a 
possibly nonlinear operator which is considered as a small perturbation. 
Let 



[/"(,(*!, i 2 ) = Texp / A(t)dt 
Ju 



be the evolution operator of the non-perturbed equation. Let us find 
the series for the evolution operator U (ti, t 2 ) of the perturbed equation 
(87) from time t\ to time t 2 by powers of the perturbation B. To this 
end, let us use the formula 

(89) U(t l7 t 2 ) = Uoih, t 2 ) + / 2 U (t, t 2 )B(t)U(t 1 ,t) dt. 

Jti 

Iterating this formula, we shall find a decomposition of the operator U, 
in the general case, as a sum over trees, on whose vertices the terms 
of the Taylor series of the operator B(t) stand, and on the edges the 
operators U stand. In the particular case when the operator B(t) is 
linear we obtain the formula 

(90) U(t!,t 2 ) = U (t ,t 2 ) (Texp £ 2 B(t)d?j U (t , h)' 1 , 
where 

(91) B(t) = U (t ,t)- 1 B(t)U (t ,t). 

We have already used this formula in the derivation of the relation (73). 
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2.4.2. Formal perturbation series for the Heisenberg equation in the tp A 
model. Let us apply this theory for the Heisenberg equation (58) in the 
ip A model. Denote the evolution operator of the Heisenberg equation 
for free field from time t to time t\ by Vo(£o,^i), so that formally we 
have 

(92) V (to, ti)$ = U (t , h) * $ * U (to, h)-\ 
where 

(93) U (t ,t 1 )=exp t -^^-H . 

(Recall that this formal expression does not exist in the Weyl alge- 
bra.) Then, by (90), the perturbation series for the evolution operator 
V(t 1 ,t 2 ) of the Heisenberg equation in the (p 4 model is given by the 
formula 

(94) V(h, t 2 )$ = V (t , t 2 ) [P(h,t 2 ) * V (t , h)- 1 ® * P^, t 2 )-\ 
where 

P{tiM) = UoiUMY 1 * U{hM) * £/o(*o,ti) 

(95) r t2 r 1 

= Texp^ J — gu(t,x) 4 /4\dtdx. 

The coefficient before g N of the latter series equals 

(96) J ^ h \N 4[ N m Tu (Hi)^m) 4 * ■ ■ ■ *u(t (N) ,x {N) ) 4 Y[dt ii) dx {ih 

where integration goes over the strip t\ < < t 2 . Absolutely the same 
integral describes the perturbation series for the evolution operator 
between any two space-like surfaces, but the integration goes over the 
strip between these surfaces. 

2.4.3. Feynman diagrams. Let us compute the expression under the 
integral in the Weyl algebra. To this end, one must firstly find the 
formula for the product of N elements of the Weyl algebra. This is left 
to the reader, starting from the case N = 3. Let us give the answer 
for the expression under the integral (96). It equals the sum over the 
Feynman diagrams, i. e. over the 4-valent non-oriented graphs with N 
vertices, and to each graph one assigns an element of the Weyl algebra 
according to the following rules: 

1) to each vertex one assigns the factor (ih)" 1 ; 

2) to each external tail (i. e. to an edge with one vertex) one assigns 
the factor u(t^, X(j)), where i is the number of the vertex; 
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3) to each edge with two vertices % and j one assigns the factor 

ih 

(97) - —T{u(t (i) , x (i) ), u(t U) , x (i) )} = -ihD(x {i) - x U) ), 

where D{x) is certain Green function of the Klein-Gordon equation, 
whose Fourier transform equals 

(98) D(p) = VV^-^ 

(W is the Cauchy principal value); 

4) to the whole diagram one assigns the factor 1/M, where M is the 
number of symmetries of the diagram, i. e. permutations of the vertices 
and the edges of the diagram preserving the graph. 

After that all factors are multiplied (in the usual sense, and not in 
the sense of *-product). 

Thus, the power of the number h for a Feynman diagram equals 
the difference between the number of internal edges and the number 
of vertices, i. e. it equals to the number of independent loops in the 
diagram minus the number of its connected components. 

We see that, for example, in the case of multiple edges the expression 
under the integral contains the square of the function D(x). This is 
a distribution with singularities on the light cone, and its square is 
non-integrable, say, for n — 3, because integral of the square of the 
expression (98) diverges at large momenta. Hence the perturbation 
series is given, in general, by divergent integrals. 

But in the tree approximation (sum over diagrams without loops) 
we formally obtain, from the perturbation series for the Heisenberg 
equation, the perturbation series for the evolution operator of the non- 
linear classical field equation 

(99) Uu -m 2 u = gu 3 /3\. 

The check of this statement is left to the reader as a useful exercise in 
perturbation theory. 

2.4.4. An attempt to define dynamical evolution in quantum field the- 
ory. The next attempt to "quantize fields" could be an attempt to con- 
struct the dynamical evolution in quantum field theory using, for each 
space-like surface, some non-commutative deformation of the algebra 
of functionals on the phase space with the Poisson bracket. Similarly 
to the linear case, in which we have chosen the deformation being the 
Weyl algebra, which admits the symplectic group of transformations, 
the required deformation in the general case could be "adapted" to 
the non-linear canonical transformations of the phase space, given by 
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the evolution operators of the Hamilton equations, i. e., of the field 
equations. To each pair of space-like surfaces one would assign an iso- 
morphism of the corresponding deformed algebras of functionals, whose 
classical limit as h — > would coincide with the isomorphism of the 
Poisson algebras of functionals, given by the transform of the classical 
evolution. 

A possible example of such deformation in the finite dimensional case 
is the so called Fedosov deformation quantization of symplectic mani- 
folds [22] . This construction uses the bundle of Weyl algebras on the 
phase space with the flat connection (Abelian connection in Fedosov's 
terminology), originating from a symplectic connection on the tangent 
bundle to the phase space. However, author's attempts to use this 
construction failed to be success. Besides that, in the finite dimen- 
sional case any Fedosov deformation is non-canonically isomorphic to 
the Weyl algebra. So it seems that for the purposes of quantum field 
theory, the Weyl algebra is the most appropriate deformation of the 
algebra of functions on the phase space, even in the non-linear case. 

Therefore the next attempt to construct quantization of fields will 
be an attempt to construct, for each pair of parameterized space- 
like surfaces Ci, C 2 , an isomorphism of the corresponding Weyl al- 
gebras Wc x — > Wc 2 - If Ci, C2 are two parameterizations of one and 
the same space-like surface, then this isomorphism should coincide 
with the action of the change of variables s on functions w l (s), p l (s). 
(Here u % (s) are transformed like functions, and p l {s) like densities.) 
For three space-like surfaces Ci, C 2 , C 3 the isomorphism Wq x — > Wc 3 
should coincide with the composition of isomorphisms — > Wq 2 and 
Wc 2 — > Wc 3 - The family of isomorphisms should be also symmetric 
with respect to the symmetry group of the theory (in the case of ip 4 
model this is the Poincare group, i. e. the group containing the Lorentz 
transformations and the parallel translations). Finally, the classical 
limit of the isomorphism Wc x — > Wc 2 as h — > should coincide with 
the isomorphism of Poisson algebras of functions given by the classical 
evolution. 



2.4.5. Dynamical evolution and perturbation theory. The subtraction 
program. Let us try to construct the dynamical evolution, as described 
in the previous subsection, for the model <p 4 , n = 3, in the framework 
of perturbation theory. Here the main idea should be the physical 
idea, due to Bethe, exposed at the beginning of the Introduction to 
Bogolyubov-Shirkov's book [13], the idea which lead to the renormal- 
ization program. Let us recall this idea in our context. Assume that 
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the required dynamical evolution exists and describes real physical pro- 
cesses for interacting fields. But in the framework of perturbation the- 
ory, we can obtain only approximations of some order with respect to 
the coupling constant, which by themselves can give divergent quan- 
tities, because the field by itself, without interaction, has no physi- 
cal sense. And the quantities which do have physical sense, such as 
the assumed dynamical evolution, can be given in perturbation theory 
by divergent expressions. For example, this means that the quantum 
Hamiltonian of the "right" dynamical evolution equals the classical 
Hamiltonian plus corrections in perturbation theory, which can be in- 
finite. The formal purpose, however, is to construct with the help of 
these heuristic constructions a "real" family of isomorphisms of Weyl 
algebras, as pointed out in the preceding subsection. 

Thus, the main idea will be an attempt to "subtract infinities from 
the perturbation series", so as to obtain convergent integrals and so 
that this subtraction of infinities have the heuristic sense of adding in- 
finite summands to the quantum Hamiltonian, which yields a family 
of (finite) isomorphisms of the Weyl algebras. Note that if we restrict 
ourselves by the space-like surfaces t = const, then we just look for a 
one-parametric group of automorphisms of the Weyl algebra of func- 
tional <&(u*(x), p l (x)). 

2.4.6. Diagram rules in the p '-representation. In order to subtract in- 
finities from the integrals corresponding to Feynman diagrams, it is 
convenient first to pass to the momentum representation. Let us state 
the rules of writing integrals in the p-representation. 

1) To each internal edge one assigns some orientation and a 4- momentum 
p, after which one assigns the factor —ihD(p) (98). 

2) To each external edge one assigns the orientation from the vertex 
outside and a 4-momentum p, after which one assigns the factor u{p). 

3) To each vertex one assigns the factor 

(100) (i/l) _1 x(±P(l) ± P(2) ± P(3) ± P(4)), 

where is the momentum outgoing from the vertex along the i-th 

edge (the sign plus is taken if the edge is oriented outside of the vertex, 
and the sign minus in the opposite case); x{p) is t ne Fourier transform 
of the characteristic function x( x ) °f the strip between the space-like 
surfaces, i. e. x( x ) = 1 if a: belongs to the strip and x( x ) = otherwise. 

4) To the whole diagram one assigns the symmetry factor 1/M. 
After that all the factors are multiplied, and one integrates over all 

the momenta p. 
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2.4.7. The "fish" diagram. Let us first consider the simplest one-loop 
diagram "fish" with two vertices (Fig. 1). To it the following integral 




Figure 1. 

corresponds: 

hrvn / xim +P( 2 ) + fc)x(P( 3 ) +P( 4 ) - fc) A, w ,„ 

(101) J (j? - m*)((k - p)* - m*) ll«(P(0)*(0*d*. 

in which we have omitted for shortness the constant and the signs W. 
This integral logarithmically diverges for large p. The divergence for 
large k is absent because of the oscillating behavior of the numerator of 
the fraction. The divergence with respect to pu\ is also absent because 
the distribution u(p(i)) is supported on the mass surface p 2 ^ = m 2 and 
rapidly decreases at infinity, since the function u(t, x) rapidly decreases 
at infinity in space directions. 

Note that the function x(p) satisfies the identity 

(102) J x(q(i) + k)x(q(2) ~ k)dk = x(q {1) + q {2) ), 

which is obtained by Fourier transform from the equality y 2 — X- If 
we subtract from the fraction under the integral the fraction 

( 1Q3 ) X(P(1) + P(2) + k)x(p(3) + P{A) - k) 

(p 2 — m 2 ) 2 

then the integral becomes convergent. Heuristically, from the initial 
integral we thus subtract the infinite expression 

/dp f 
(p2 _ P m 2)2 J +P(2) +P(3) +P(4))Y[u(P(i))dp {i) , 

which corresponds to subtraction from the Hamiltonian of the infinite 
term 

., 9 f dp u 4 
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Hence at this level the subtraction program gives a correctly defined 
family of isomorphisms of the Weyl algebras, satisfying all the necessary 
requirements. This family of isomorphisms, however, is defined not 
uniquely, but only up to adding a finite summand ihg 2 cu A /'i\ to the 
Hamiltonian. 



2.4.8. The two-loop diagram. Let us now consider the two-loop dia- 
gram with two vertices (Fig. 2). To it the following integral corre- 




FlGURE 2. 

sponds: 

(106) I X(P(D + fc)x(P(2) ~ k) 

J (p 2 — m 2 ) (q 2 — m 2 ) ((k — p — q) 2 — m 2 ) ' 

in which we have omitted for shortness the constant before the inte- 
gral, the symbols VV and the factor flLi ^{P{i))dP{i)- The integral is 
divergent for large p, q. To make it convergent, one can, for instance, 
subtract from the fraction 

(107) ' 



(p 2 — m 2 )(q 2 — m 2 )((k — p — q) 2 — m 



2\ 



its Taylor polynomial with respect to k at k = of the second order, 
i. e. the terms of the zeroth, first and second order in the Taylor devel- 
opment. Then the remainder will be an integral of partial derivatives 
with respect to k of the third order, and it is not difficult to see that it 
would give the convergent integral instead of (106). But in this process 
in the numerator the integrals 



(108) 



j hx(P(i) + k)x(p(2) - k)dk, 
J kikjx(p{i) + k)x(p(2) - k)dk, 
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will occur, which are the Fourier transforms with respect to the variable 



expressions are not defined as distributions. Here the problem is that 
the characteristic function x( x ) i s n °t differentiate. 

Thus, we see that our program of defining dynamical evolution in 
quantum field theory fails on the two-loop diagram. 

2.4.9. Dynamical evolution in the quasiclassical approximation. How- 
ever, in the one- loop approximation, i. e., in the sum over the diagrams 
with no more than one loop, the program of defining dynamical evolu- 
tion works well. We come to the following theorem. 

Theorem. In the y? 4 model of quantum field theory in four di- 
mensional space-time the dynamical evolution exists in the one-loop 
approximation of perturbation theory. 

This theorem means that to each pair of space-like surfaces Ci, C 2 
one can assign, with the help of the subtraction procedure, an element 
of the Weyl algebra W of the type 



up to 0(h), the evolution operator of the classical field equation (99), 
and 



Besides that, as the surface C\ tends to t = — 00 and the surface C 2 to 
t — 00, the element Pc x ,c 2 tends to the S-matrix up to o{h) (see the 
next Subsection). 

Proof of this theorem is based on the fact that the only one-loop 
diagrams giving divergent integrals are the diagrams containing the 
diagram "fish" from 2.4.7. Making the same subtraction procedure 
with them as with the "fish" diagram, we shall obtain the required 
element Pc lt c 2 - 

This theorem is in accordance with the results from the book [2] by 
Maslov and Shvedov, who constructed complex germ in quantum field 
theory using the Bogolyubov S'-matrix (regarding this S'-matrix see the 
next Subsection). 

2.4.10. The scattering matrix. Thus, for the two-loop diagram the sub- 
traction program meets the difficulty that the characteristic function 
x(x) of the strip is not differentiable. Let us slightly change the view- 
point on the dynamical evolution, and let us look not for a family of 
isomorphisms of the Weyl algebras of space-like surfaces, related with 





(110) 



Pdfia = Pc 2 ,C 3 * PC U C 2 + 0(h). 
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the integral (95), but for one element of the Weyl algebra, playing 
the role of the evolution in the whole space-time and related with the 
integral 

(111) Texp J J ^(t,x)w(t,x) 4 /4!rftrfx, 

where g(t,x) is a smooth function, say, with compact support. In 
other words, let us change the function gx( x ) m our considerations by 
a differentiate function g(x), and consider the Lagrangian 

^/ „ 1 m 2 9 q(x) a 

(112) F(x tl , u, u x ?) = -u x ,u x ^ - —u 2 - ■ 

The perturbation series for this Lagrangian is given by the integral 
(111). 

If we develop this integral according to the rules from 2.4.3 and 2.4.6, 
i. e. develop the summand 
(113) 



/ 



(ih)N4\N m 9{x(i)) ■ ■ ■ g(x (N) )Tu(x (1) Y * ... * u(x (N) ) 4 Yl dx (l) , 



where integration goes over the whole space-time, then the obtained 
integrals, related with the Feynman diagrams, will exactly coincide 
with the Feynman integrals from Bogolyubov-Shirkov's book [13], with 
the only difference: instead of function (98), in the Feynman integrals 
the propagator 

(H4) D c (p) = L— _ 

p 2 — m 2 + ie 

stands. This propagator differs from D(p) by a multiple of the delta- 
function S(p 2 — m 2 ), hence this difference does not affect on the di- 
vergences at large momenta. Hence we can apply to our integrals 
the subtraction procedure from the Bogolyubov-Shirkov's book (the 
Bogolyubov-Parasyuk theorem). In fact, we have already begun to 
apply it in 2.4.7, 2.4.8. 

Having applied this procedure, we will obtain an element P(g) of 
the Weyl algebra Wq, which is a formal series over the powers of the 
function g(x) and which is defined not uniquely, but only up to adding 
finite terms to the Lagrangian. Conjugation by the element P(g) in 
the Weyl algebra gives, up to 0(h), the perturbation series for the 
evolution operator of the classical field equation 

(115) Du(x) - m 2 u(x) = g(x)u 3 (x)/3\ 

from t = — oo to t = oo. 
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Further, consider the operator in the Fock space corresponding to 
the element P(g) (see 2.3.3). Denote it by S{g). We state that the 
operator S(g) is exactly the S-matrix constructed in the book [13] by 
Bogolyubov and Shirkov. This S-matrix is obtained from the element 
P(g) by the change everywhere of function D(p) by the Feynman prop- 
agator D c (p), the *-product of functionals by the composition of oper- 
ators, and the usual product of functionals (for example, m 4 (x)) by the 
normally ordered product of operators. Indeed, by the Wick theorems 
(see [13]), operations with products of functionals in the Weyl algebra, 
such as *-products with the sign T, exactly correspond to operations 
with normally ordered products of operators in the Fock space, with 
the only difference: the function D(x) should be replaced in these for- 
mulas by the function D c (x). And the subtraction procedure in the 
Weyl algebra exactly goes to the subtraction procedure for operators 
in the Fock space. 

Thus, the operator S(g) is the Bogolyubov S -matrix, or the scattering 
matrix. This S'-matrix satisfies the Lorentz invariance, unitarity and 
causality conditions. In the Weyl algebra these conditions go to the 
corresponding conditions for the element P(g)- The Lorentz invariance 
condition is obvious: 



for a Lorentz transformation L. The unitarity condition means that 



Finally, the causality condition states that for two functions g±(x) and 
g 2 {x), coinciding for t < t , the element P(gi)*P(g 2 )~ 1 does not depend 
on the behavior of the functions g±, g 2 for t < t . These conditions are 
the natural substitutes of the conditions for the dynamical evolution. 
In particular, the causality condition is the natural substitute of the 
condition of dependence of the evolution operator of a linear differential 
equation on the coefficient functions of this equation. In Bogolyubov- 
Shirkov's book [13] it is shown that these conditions define the S'- 
matrix uniquely, up to adding finite terms to the Lagrangian. Hence 
it is natural to postulate the existence of the elements P(g) and S(g) 
outside of the framework of perturbation theory. 

The physical sense of the S- matrix is the following. Assume that, 
as the function g(x), remaining a function with increasing compact 
support, tends to the constant function g = const (the adiabatic inter- 
action switch off), the elements S(g) and P(g) tend to some elements S 
(the physical S-matrix) and P. Then the square of the absolute value 
of the matrix element (86) (with $ = P) of the physical S-matrix is 



(116) 




(117) 



P{g)*P(g) = l. 
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the density of probability of the event that iV colliding particles, flying 
before collision with 4- momenta . . ., P(jv), after collision turn into 
N' particles flying away with 4-momenta p',^, . . ., p'/ N ,y 

Up to o(h) the element P gives the operator Pc ll c 2 °f quasiclassi- 
cal dynamical evolution from C\. t = — oo to C^- t — oo, see 2.4.9. 
The definition of dynamical evolution (up to o(/i)), suitable outside 
perturbation theory, can be found in 2.4.4. 

Let us also comment on the non-uniqueness of S-matrix. Actually 
it depends not only on the initial Lagrangian, but also on effective 
parameters, for example, on effective mass and effective coupling con- 
stant, which are computed from the S'-matrix and which already define 
it uniquely. And any change of effective parameters is equivalent to 
certain change of parameters of initial Lagrangian. 

In a similar manner one constructs the apparatus of the operator 
and scalar Green functions, with the help of the Lagrangian 

(118) F(x M , u, u x ») = -«x"M % - — u 2 —u 4 -j(x)u. 

Their construction, including the subtraction procedure, does not yield 
new difficulties, and is made similarly to what is done in the Bogolyubov- 
Shirkov's book. 

Note that the above constructed apparatus of the S'-matrix and con- 
ditions on it are analogous to the scattering theory in the theory of 
partial differential equations, where, given the coefficient functions of 
the equation, one is required to determine the properties of the evolu- 
tion operator from t = — oo to t = oo. 
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